Nonextensive scaling in a long-range Hamiltonian system 
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The nonextensivity of a classical long-range Hamiltonian system is discussed. The system is 
the so-called ct-XY model, a lattice of inertial rotators with an adjustable parameter a controlling 
the range of the interactions. This model has been explored in detail over the last years. For 
sufficiently long-range interactions, namely a < d, where d is the lattice dimension, it was shown to 
be nonextensive and to exhibit a second order phase transition. However, conclusions in apparent 
contradiction with the findings above have also been drawn. This picture reveals the fact that there 
are aspects of the model that remain poorly understood. Here we perform a thorough analysis, 
essaying an explanation for the origin of the apparent discrepancies. 
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I. INTRODUCTION 



Systems of many particles interacting via long-range 
forces, although ubiquitous, are not fully understood (see 
for instance [1]). Special interest in such systems has 
arisen recently in connection with the extension of stan- 
dard statistical mechanics proposed by Tsallis [2]. As 
a prototype to study the dynamics and thermodynam- 
ics of long-range systems, both in equilibrium and non- 
equilibrium situations, a dynamical model with an ad- 
justable interaction range has been introduced [3]. The 
model consists in N interacting rotators moving on par- 
allel planes and located on a periodical d-dimensional 
hypercubic lattice with unitary spacing. Each rotator is 
fully described by an angle < 0$ < 2ir and its conjugate 
momentum Lj. The dynamics of the system is ruled by 
the Hamiltonian 
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K + V , 
(1) 



where the coupling constant is J > (we restrict our 
study to the ferromagnetic case where interactions are 
attractive) and, without loss of generality, moments of 
inertia equal to one are chosen for all the particles. Here 
Tij measures the minimal distance between rotators lo- 
cated at the lattice sites i and j. One can associate to 
each rotator a "spin vector" m.j = (cos6>j,sin#i), which 

allows to define an order parameter m = -^fX^=i m i- 
The Hamiltonian (1), describing a classical inertial XY 
ferromagnet, is usually referred to as a-XY model. It 
includes as particular cases the first-neighbor (a — > oo) 
and the mean-field (a — 0) models. Note that this is 
an inertial generalization of the well known XY model 



of the statistical physics of magnetism: the time evolu- 
tion is given by the natural dynamics governed by the 
Hamilton equations. 

This prototype of complex long-range behavior has 
been thoroughly explored in the last few years (see for in- 
stance [3-5]). It has been shown that the model presents 
nonextensive behavior for a < d [3]. In that domain of 
a it displays a second order phase transition. This result 
has been exhibited first by means of numerical compu- 
tations for the one-dimensional (ID) case [4] and later 
through analytical calculations for arbitrary d using a 
scaled version of the Hamiltonian H [5]. However, a re- 
cent work [6] draws conclusions that are in disagreement 
with the previous findings, claiming that the model is 
extensive for all a and that there is no phase transition. 
This apparent contradiction helps to put into evidence 
that there are aspects of the model that remain obscure. 
The lack of a comparitive study as well as of a discussion 
on the origin of the discrepancies motivates the present 
work. It is the purpose of this paper to review and com- 
plement previous results to elucidate the question. 

In order to do that, we use the following methodology. 
We start by solving the equations of motion associated 
to Hamiltonian H: 



L, 



dH 

dH 

d6i 



= U 



= -J 



(2) 



y-^ sm( 



i=l,...,N. 



Numerical integration is performed by means of a sym- 
plectic fourth order algorithm [7] using a small time 
step to warrant energy conservation with a relative error 
smaller than 10~ 5 . Equilibrium properties are analyzed 
by means of time averages (computed after a transient) 
that allow to mimic mirocanonical averages. In Ref. [6], 
numerical results for the canonical ensemble were ob- 



*e-mail: celia@cbpf.br 



1 



taincd through standard Monte Carlo simulations. Due 
to ensemble equivalence [8], both methods are expected 
to yield the same macroscopic averages at thermal equi- 
librium. Simulations will be supplemented by analytical 
considerations. 



II. EQUILIBRIUM THERMODYNAMICS OF 
THE a-XY MODEL 

Along this paper we will consider Hamiltonian (1) al- 
though many related works in the literature refer to a 
modified version of this Hamiltonian such that the in- 
teractions are scaled. Since there is a correspondence be- 
tween both descriptions, we will discuss this point to take 
profit of all the pertinent results in the literature. To con- 
struct the scaled Hamiltonian, let us call it H, the cou- 
pling coefficient J in H is substituted by J/N, where N is 
the upper bound of the potential energy per particle, that 
depends on N, a, and d according to N = j^y ^ 5^7^- 

i j^i 

In the large N limit one has [4,9] 

(- N 1 - a / d 0<a<d 
N(N,a/d)~l In A a = d (3) 

I Q(a/d) a>d 

with 9 a function of the ratio a/d only, that goes to 2 as 
a/d goes to infinity. For a < d, N depends strongly on 
A. Then the representation given by H may be consid- 
ered artificial, since the microscopic coupling coefficient 
becomes A-dependcnt, that is, becomes fed with macro- 
scopic information. Anyway, the thermodynamics and 
the underlying dynamics of H can be trivially mapped 
onto those of H by transforming energy-like quantities 
through E <-» E/N and characteristic times (as long as 
moments of inertia remain unitary) through r <-» rA^ 
[3]. The usual preference for the scaled form H comes 
from the fact that the thermodynamic limit of E/N is 
always finite and no further scalings of either thermody- 
namical or dynamical quantities are needed. 

It has been analytically shown [5], through canoni- 
cal calculations performed with Hamiltonian H , that the 
thermodynamics of systems with a < d, at the final ther- 
mal equilibrium, is equivalent to that of its mean-field 
version (the so-called Hamiltonian Mean Field (HMF) 
[10]). Such systems display a second order phase tran- 
sition, from a low-energy ferromagnetic state to a high- 
energy paramagnetic one, at a certain critical energy per 
particle e c — U c /N — 0.75 J. This important result for H 
analytically confirms the previous findings [4] for the ID 
case of the original Hamiltonian H, just by taking into 
account the simple mapping between H and H. Since the 
equilibrium results of the HMF are universal for a < d, 
we will summarize them. In terms of the magnetization, 
Hamiltonian H leads to the following caloric curve 



A JA N 2 , . , hWJm) 

U = -= + l-m 2 , with to — — — ^ '-, 

2/3 2 L J Io(PJm) 

(4) 

where f3 = 1/f (being f = 2(A) /A the temperature 
and having set the Boltzmann constant ks = 1) and /„ 
are the modified Bessel functions of order n. The consis- 
tency equation from which the magnetization is extracted 
can be found for instance through canonical calculations 
[10]. It has a stable solution to = for (3 J < 2 while, 
for f3J > 2, the zero magnetization solution becomes un- 
stable and a non-vanishing /3-dependent stable solution 
arises. From (4), it is clear that the critical value j3 c J = 2 
corresponds to e c = 0.75J. Notice in Eq. (4) that, as 
to 2 < 1 and the inverse temperature (3 does not depend 
on A, then the large A limit of E/N is always finite. 

We will analyze the size dependence of thermal aver- 
ages. We will focus on the range 0<a<loflD lattices 
governed by Hamiltonian H. In Fig. 1(a), the average 
magnetization per particle (to) is represented as a func- 
tion of the energy per particle U/N, for a — 0.5 and 
different system sizes. Clearly, the energy per particle 
at which the system becomes disordered, i.e, at which 
the magnetization vanishes up to finite size corrections, 
grows with the system size. In Fig. 1(b), the same data 
are represented as a function of U/NN. Through this 
scaling, all data sets tend to the same curve in the ther- 
modynamic limit, as it has already been shown previ- 
ously [4]. In the inset of Fig. 1(b), a plot (to) vs. A, 
for U/NN = 1.4 and two values of a, illustrate that 
the magnetization in the high energy regime decays with 
the system size as 1/VN. Additionally, data sets for 
A = 128 and different values of a e [0, 1) were included 
in Fig. 1(b) to show that the curve of magnetization vs. 
U/NN is the same for any a-XY system with < a < 1, 
up to corrections of order l/y/~N. In particular, the uni- 
versal curve coincides with the one for the HMF model 
(a = 0), given by Eq. (4), once taken into account the 
mapping E <-> E/N. Everything in agreement with the 
analytical results of Ref. [5]. However, concerning the 
phase transition, there is a risk to fall into an endless 
rhetorical discussion. Strictly speaking, there is no ferro- 
magnetic transition, because the critical energy per par- 
ticle U c /N — 0.75 J N is divergent, as asserted in [6]. 
Nevertheless, the limit A — > oo, despite being an ideal- 
ized situation, must reflect the behavior of finite but large 
systems in order to be meaningful. Ultimately, we are in- 
terested in finite-size systems, as real systems are. For 
finite-size a-XY systems, with a < 1, as those that were 
simulated in this work, one can distinguish two regimes: 
One, at low energies, where the system is ordered with a 
magnetization significantly different from zero and inde- 
pendent from the system size, and another, a disordered 
one, with magnetization of order l/\/A. A good, repre- 
sentative, thermodynamic limit, reflecting this situation, 
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can be denned by means of an appropriate scaling, the 
one allowing data collapse. In that case, it results a fi- 
nite critical energy, U C /NN = 0.75J, which plays the 
same role as the critical energy per particle in an exten- 
sive system. 
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FIG. 1. (a)-(b) Average magnetization per particle (m) as 
a function of the total energy U for different sets of param- 
eters (a,N). The inset in (b) presents the magnetization as 
a function of 1/VN for U/NN = 1.4. (c)-(f) Caloric curves. 
In all cases full lines correspond to the mean-field analytical 
result given by Eq. (4) and, dotted lines are guides to the 
eyes. Dashed lines correspond to U/N = 1//3 in (c) and to 
a straight line of slope 1/2 in (d). Symbols (defined as in 
(b)) correspond to averages over 10 samples, computed over 
a time interval of order 10 3 , after a transient (t w 10 3 ) has 
elapsed and starting from "water-bag" [10] initial conditions. 
We have set J = 1. The lattice dimension is d — 1. 

Let us analyze the dependence of the mean kinetic en- 
ergy on the total energy for different system sizes and 
different values of a. In Fig. 1(c) we represent the 
data as done in Ref. [6], that is U/N vs. 0, where 
1//3 = T = 2{K)/N. In agreement with Ref. [6], perfect 
data collapse occurs for total energies scaled with the sys- 
tem size N, as can also be appreciated in the alternative 
representation of the same data exhibited in Fig. 1(d). 
However, these plots are restricted to very low energies. 



If one extends the range of energies plotted (Fig. 1(e)), 
it becomes clear that data collapse does not hold any 
more through the TV-scaling. Whereas, as before, it is 
the TV TV-scaling the one which leads to data collapse in 
the full energy range (Fig. 1(f)). As an aside comment, 
note that, because the relation U ~ 2 (if) holds at low 
energies, data collapse would occur in that regime for any 
arbitrary scaling by N 1 , with 7 6 3?. In particular, this 
is true for 7 = 1, as plotted in Fig. 1(d) (hence Fig. 1(c) 
at low energies) and for 7 = 2 — a, as in Fig. 1(f) at low 
energies. 

One can understand what is going on as follows. For 
very low energies, the dynamics is dominated by the 
quadratic terms of the potential. Thus, the system can 
be seen as a set of almost uncoupled harmonic oscilla- 
tors (normal modes). One can also think of particles in 
a mean-field, a description that is exact in the infinite- 
range case. The particles effectively interact not through 
the full mean-field m but only through its fluctuations. 
If the mean-field were constant it would play the role of 
an external field and there would be no interactions. At 
low energies, where m is almost constant, the residual 
or effective interaction, that is the component coming 
from the fluctuations of m, is small. This is consistent 
with the normal modes view, where interactions are very 
weak too. Therefore, in the limit of very low energies (as 
well as in the limit of very high energies) the system be- 
comes non- interacting (hence, integrable). While at high 
energies, i.e., above the critical value, one has almost 
non-interacting rotators; at low energies, i.e., close to 
the ground state, one has almost non-interacting normal 
modes. Then, at low energies, from the virial theorem, 
the result (K) ~ (V") arises trivially. The consequent re- 
lation U ~ 2 (if) = N/(3 indicates that the energy is ex- 
tensive. A natural result since the interaction terms are 
not strong, contrarily to what was asserted in [6]. How- 
ever, as the energy increases and anharmonicitics grow, 
the correct scaling choice is no more that of an extensive 
system, as becomes evident in Fig. 1(e). Data collapse is 
actually obtained through the scaling by TV TV, as shown 
in Figs. 1(b) and 1(f). Moreover, this data collapse is 
expected to be universal for any a £ [0,d) [5]. Hence, at 
criticality, we have the nonextensive behavior U c oc JNN 
and also 1/(3 C = T c oc JN. 

Let us review the whole picture from the viewpoint 
of canonical ensemble calculations. We will consider 
the case a = 0, but although tricky, a generalization 
to arbitrary a £ [0, d) could be analytically performed 
[5]. The partition function of Hamiltonian (1) when 
a = is given by the following integral over phase space 
Z = J n*Li dljdOj exp(— (3H) = ZnZy, which factorizes 
into the kinetic and potential contributions 
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where Zy has already been transformed by means of 
the Hubbard-Stratonovich trick and G{y) — — -^lny + 

2 

2 p JN — l n ^o(y)- The derivation is the same followed 

in [10] for the scaled Hamiltonian H , apart from an N- 
scaling that does not affect the procedure. The integra- 
tion can be performed by means of the Gaussian approx- 
imation around the point y verifying G'(y ) = 0, that 
is, y Q ~ (3JNh(y o )/I {y o ) and G"(y ) > 0. In our case, 
the total energy results 

d\nZ N JN 2 2 , 
U = -^T=2l3 + —( l - m) > (6) 

with 

hipJNm) 

Tfl == 

hifiJNmY 

in correspondence with Eq. (4). For large (3JN, from 
the consistency equation, one has m 2 ~ 1 — -^j-^ for the 
stable solution, an approximation that is equivalent to 
considering y Q ~ (3JN — 1/2, as done in Ref. [6]. In fact, 
substitution of the above approximate expression for m 2 
in (6), gives U ~ N//3. Again one obtains that at low 
temperatures the energy is extensive. However, the ap- 
proximation above is no longer valid as (3 decreases. In 
this case, long-range couplings become effective and the 
nontrivial nonextensive behavior comes out. Then, the 
energy no longer scales with N and one has to consider 
the more general Eq. (6) . An analysis as that performed 
in Ref. [6], restricted to the very low temperature regime, 
misses most of the rich physics of the long-range interact- 
ing rotators. Of course, this discussion is meaningful as 
soon as N is not excessively large. Recall that 1 j(i c <~ JN 
for generic a, hence 1//3 C ~ JN for a = 0. Then N has 
to be large enough so that the thermodynamic limit is 
a reasonable approximation but not so large as to drive 
the temperature scale out of a realistic range. 

III. FINAL REMARKS 

A double sum as in (1) indicates that, for interaction 
ranges < a/d < 1, the total energy U may grow as 
A 7 , with 7 > 1, as occurs in the regimes of the a — XY 
model where long-range couplings become relevant (see 
also [11]). Therefore, the large N limit of U /N is not well 
defined, in fact, the energy per particle diverges when 
N — > oo. In that case, the energy is a nonextensive 
quantity [12]. Many systems in nature also display such 
kind of behavior, as illustrated by Thirring in the context 
of a discussion on the stability of matter [13]. In those 
cases, it is sometimes said that the thermodynamic limit 
does not exist. However, a proper thermodynamic limit 



can be effectively achieved by introducing a suitable N- 
dependent factor N* <~ TV 7-1 such that the large TV limit 
of U/NN* results well defined [14]. Concerning critical- 
ly, for the a — XY model, in the thermodynamic limit, 
there is no phase transition in the sense that a transition 
does never occur at a finite energy per particle. However 
for finite-size a-XY systems, with a <G [0, d), one can 
distinguish two regimes: An ordered one at low ener- 
gies and a disordered one above a "critical" energy that 
increases nonextensively with the system size (see Fig. 
1(a)). Then, a different limit appears to be the relevant 
one. Indeed, application of an appropriate rcgularization 
procedure, namely, further scaling by N* — N, allows to 
display a transition. By means of that scaling, a finite 
critical energy, U C /NN = 0.75J, can be defined. In this 
way, a thermodynamic limit, representative of the be- 
havior observed for large N (although not exceedingly), 
which is not limited to the low energy clustered regime, 
is obtained. 
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